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1. Introduction 

The production of hadronic jets in proton-(anti)proton collisions is one of the most basic pro- 
cesses at hadron colliders. Key processes include those with a W, Z, or Higgs boson together with 
jets. Many models for new short distance physics produce colored particles, which would exper- 
imentally be observed as jets of energetic hadrons. For example, in supersymmetric models the 
squarks or gauginos typically decay through a chain of (colored) particles into the lightest super- 
symmetric particle. Another source of jets are the initial-state quarks and gluons taken out of the 
incoming protons, which radiate before entering the hard collision. 

The underlying hard process we are ultimately interested in is always accompanied by QCD 
corrections contributing at the various scales that are present in the process: virtual corrections to 
the hard process itself, real and virtual corrections from initial-state and final-state radiation, etc. 
The QCD radiative corrections to processes with jets are typically large. One reason is that jet 
selection cuts can be sensitive to additional soft and collinear emissions. In perturbation theory 
this sensitivity to lower scales manifests itself via large logarithms of the form a" ln m \ljQ with 
m < 2n. Here, Q is a large scale of the order of the partonic center-of-mass energy and }X is a 
low scale associated with the definition of the fmal-state jets of the order a few tens of GeV. In 
the limit jJ. <^ Q these large logarithms degrade the fixed-order perturbation series, resulting in 
large perturbative uncertainties. This can be avoided (or at least alleviated) by resumming the 
leading (and subleading) towers of logarithms to all orders in a s . In practice, this results in smaller 
perturbative uncertainties, and hence more reliable predictions for the cross section in question. 

In inclusive jet measurements, one requires a minimum number of hard jets and sums inclu- 
sively over additional emissions. In this case, one is less sensitive to lower jet scales. On the other 
hand, for exclusive jet cross sections, where one requires a certain fixed number of jets in the fi- 
nal state, one is explicitly sensitive to the jet resolution scale through the veto on additional jets. 
Such exclusive jet measurements play an important role in Higgs measurements and new-physics 
searches at the LHC. 

An effective-theory framework to disentangle the relevant scales in jet production and resum 
the associated logarithms of ratios of these scales via renormalization group evolution (RGE) is 
provided by Soft-Collinear Effective Theory (SCET) [Jl], g, g, g, g, §]. In SCET the soft-collinear 
limit of QCD is manifestly implemented at a Lagrangian and operator level using a systematic 
power expansion. The QCD corrections appearing at the hard-interaction scale, i.e., away from any 
infrared singular limits, which contain the process-specific details, are incorporated via matching 
from QCD onto SCET. On the other hand, the infrared-sensitive QCD dynamics below the hard- 
interaction scale that describes the collinear radiation within the jets and soft interactions between 
jets is contained in the effective theory. 

2. Resummation of Exclusive Jet Cross Sections in SCET 

Consider a process with N final-state jets and L leptons, photons, or other non-strongly inter- 
acting particles with underlying hard interaction 

K a {q a ) K b {q h ) -)■ K\(q\) ■■ ■ K N+L (q N+L ) , (2.1) 
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where K a p denote the colliding partons, and jq denote the outgoing quarks, gluons, and other 
particles with momenta qi. The incoming partons are along the beam directions, q^ b = x a .bPa,b, 
where x a are the momentum fractions and P a j } the (anti)proton momenta. We are interested in 



the situation where each of the partons in Eq. ( |2.1| ) produces a separate identified jet and where 
we do not allow additional hard jets coming from a hard ISR or FSR emission. In SCET, the 
corresponding exclusive jet cross section can be factorized, which leads to an expression of the 
form [|, 0,| 



da = / dx a dx b d$> N+L (q a + q b \{qi})F '({<?;}) £tr [H K (q a ,b,{qi})S^ i 



B K a B K b Y[ J 



(2.2) 



Here, A^^ + i{q a +q b , {qi}) denotes the Lorentz-invariant phase-space for the final-state particles in 



Eq. (2.1), where F({qi}) denotes the measurement made on the hard momenta of the A" signal jets 
(which in the factorization at leading order in the power expansion are approximated by the parton 



momenta qi). Equation ( |2.2| ) does not include possible contributions from Glauber exchange. 

Any dependence probing softer momenta, such as measuring jet masses or low px&, as well 
as the choice of jet algorithm or jet resolution variable, will affect the validity and precise form 



of the factorization formula Eq. ( |2.2| ). This dependence is encoded through the precise definitions 
of the A-jet soft function (describing soft radiation), the jet functions J Kj (describing energetic 
final-state radiation in a single jet) and the beam functions B Kib (describing energetic initial-state 
radiation). The parton distributions of the incoming protons, fj, are contained in the beam func- 
tions, which can be further factorized as 5; = Y,j ^ij ®fj » 10 ] . The beam, jet, and soft functions 
contain the virtual and integrated real-emission corrections in all IR-singular limits and are sepa- 
rately IR finite. 

A jet-resolution variable with particularly simple factorization properties is A-jettiness [|0|], 



which effectively provides an exclusive A-jet algorithm. The explicit form of Eq. ( |2.2| ) for the A-jet 
cross section defined using A-jettiness is known [|| 11, 12], and for A^ = the resummation has 
been carried out to NNLL [O, 14]. The resummation of exclusive jet cross sections defined using 



more traditional jet-clustering algorithms has been studied for example in Refs. [|l5|, [l6[ [l7[ [TJ, [l9|, 



20, 21]. 



The remaining ingredient in Eq. (2.2) is the hard function H K (q a ,q b , {qi}), where the sum 



over K = {x a , K b , . . . Kjv+z,} in Eq. (2.2) is over all relevant hard-interaction channels. The hard 



function encodes the dependence on the underlying hard interaction Eq. (2.1), including the hard 
virtual corrections. It is explicitly independent of the used jet definition or jet resolution variable 
and therefore does not depend on the precise form of the factorization. Once the resummation for 
a given jet observable and a given number of jets to a certain order is known, it can in principle be 
applied to any desired process. In most cases the bottleneck is then to extract the required process- 
specific information in form of the NLO hard function from existing fixed-order calculations. 



3. Matching from QCD onto SCET With Helicity Amplitudes 

Schematically, we match onto the effective Lagrangian 

-2eff = -S^SCET +2^Ck Ok , 
k 



(3.1) 



3 



Combining fixed-order helicity amplitudes with resummation using SCET 



Frank J. Tackmann 



where J;?scet is the SCET Lagrangian for soft and collinear quarks and gluons at leading or- 
der in the power expansion. The operators Ok are responsible to mediate the hard interaction in 



Eq. (|2JJ) by creating and destroying the appropriate number of external partons. The hard function 
in Eq. ( ^2| ) is given by the square of the Wilson coefficients: H = \Ck\ 2 - They are determined by 
requiring that the UV-renormalized amplitudes in the full and effective theories agree, 



*QCD — -^SCET 



£iQ<0*) SCBr . (3.2) 



The Wilson coefficients explicitly depend on the UV renormalization scheme adopted for the SCET 
operators 0%, for which we use the standard dimensional regularization together with MS. 

As in any effective-theory setup, the IR divergences of the full theory are by construction re- 
produced in the effective theory and cancel in the matching. In particular, the matching coefficients 
Ck do not dependent on the specific IR regulator used to perform the calculation. A useful choice is 
dimensional regularization for both UV and IR divergences, then all loop graphs in SCET are scale- 
less and vanish. Thus, the UV and IR divergences in SCET precisely cancel each other, and the 
bare matrix elements are given by their tree-level expressions. Including the counter term 5o(£uv) 
due to operator renormalization in the effective theory removes the UV divergences and leaves the 
IR divergences. Schematically, the SCET amplitude is thus 

^scET = iC-((0) tree + (0) lo °P) =iC[l + S (£ IR )] . (3.3) 

Since the effective-theory IR divergences, C8o{£m), nave to match those of the full theory, the 
matching coefficients in MS are directly given by the infrared-finite part of the full-theory ampli- 
tude computed in pure dimensional regularization, C = — I^qcd- 

The above simplification provided by pure dimensional regularization is well known, and was 
used for processes with multiple external partons for example in Refs. [22, 23, ^4]]. Our goal is 



to construct a general and convenient-to-use operator basis, which lets us exploit it as much as 
possible. This requires to organize the possible spin and color structures, which quickly proliferate 
when increasing the number of external legs. It should also be easy to incorporate constraints from 
charge conjugation and parity invariance and to remain fully crossing symmetric. As one might 
expect, this can be achieved by constructing the operator basis by employing the same helicity and 
color decompositions as used for the QCD amplitudes. As a result, the SCET Wilson coefficients 
will be given directly by the IR-finite parts of the QCD color-ordered helicity amplitudes. 
For each jet direction, we define two types of light-cone vectors 

nf = (l,nO, nf = (l,-ni), (3.4) 

with nj = nf = 0, • hi = 2, and = qij\q\ \ is a unit three-vector in the direction of the /th jet. At 
leading order in the power expansion a fixed-order QCD amplitude with N colored legs is matched 
onto operators in SCET with N different collinear fields, where the different collinear directions 
have to be well separated, • nj 3> A 2 for i ^ j. The SCET operators are constructed out of 
composite fields that are invariant under collinear gauge transformations |||] 



Xn.i 



(3.5) 
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Here, £,„ and A„ are «-collinear quark and gluon fields. The «-collinear covariant derivative and 
Wilson lines are denned as 



perms N " ^« 



(3.6) 



The Wilson line W„ (x) sums up arbitrary emissions of «-collinear gluons from an n-collinear quark 
or gluon, which are &(\) in the power counting. 



Using the polarization vectors defined in the standard spinor notation (see e.g. Ref. Q25P), 

we define a gluon field of definite helicity 

®f ± = -e^(n i ,h i )^ aiXi , (3.8) 

where a is an adjoint color index. For example, for = (1,0,0, 1) (and with an appropriate spinor 
phase convention) we have 

(n h ni) = ±= (0, 1 , Ti, 0) , <%f ± = -±= {X^ ± i.<i ) . (3.9) 

Similarly, we define a gg-current of definite helicity 

jt = ^(PuP ^ 1 -^^ , Pt = «k%, (3-10) 

where a, /3 are fundamental color indices, and pi denotes the large label momentum carried by 
collinear fields. The corresponding tree-level Feynman rules are 



g a ± (p)\^l\0) = 8 ab 5(p i -p), (3.11) 

m 



{q a ±(Pi)f4(Pj)\ J ff± 1°) = ^ 8aA d (Pi ~ P^ 5 (PJ ~ Pi) > ( 3 - 12 ) 



while any other combinations of helicities vanish. 

Next, we assemble the helicity fields into operators with a definite helicity structure corre- 
sponding to each external helicity configuration, 

n a\a%— a;_ia;— csv-ia^/ ~ ~ ~ ~ ~ ~ \ _ o ^2 jOi-iUi T a N -i<x N 

u ±±...(±... ; ...±) \P\iP2i--- ,Pi-l,Pi, ■ ■ ■ iPN-1,PN) — 0^ 1± ^ 2± • ■■J qi _i i± ■ " J q i}{-\,N± • 

(3.13) 

The symmetry factor S takes into account the number of identical particles. For each number of 
positive and negative helicity gluons and quark currents, there is only one independent operator. 
This is because fields of the same particle type and helicity are related by interchanging their 
momentum labels and color indices appropriately. For example, O a ^ a2 (p\ : p2) = O a _}^ [ (p2,P\) are 
not independent operators. To keep track of the minimal number of independent operators, we can 
simply order the helicity labels. 
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For a given A^-parton process, the explicit form of Eq. (3.1 ) reads 

helicity J i=\ 
configurations 

Using the above Feynman rules for the helicity fields, the resulting tree-level Feynman rule for 
jSfeg- projects out the single Wilson coefficient that corresponds to the external helicity and color 
configuration, 

<^(pi)g- 2 (p 2 )--^-H^ 

(3.15) 

Using pure dimensional regularization with MS and following the same arguments as above, we 
thus find that to all orders in perturbation theory 

C a ;r { .^) iaN (Pi,P2,---,PN-uPN) = -i^(g^ipi)g^(P2)---q- N - 1 ipN-x)q c t N iPN)) ■ (3-16) 

Next, to organize the color structures, we can decompose the Wilson coefficients into a com- 
plete basis of color-singlet structures, 

c +-(—) — 1* )h =1 Kp-i-i) 

k 

Here is a row vector of suitable color structures that provide a complete basis for all 

possible allowed color structures (but do not necessarily all have to be independent). Up to three 
colored partons, there is only a single allowed color structure 

fUfi = ; f ub =(8 ab ), f Ua P = (T^p), f ubc = (if bc ), (3.18) 

while for example for qqqq or ggqq there are already two or three independent color structures, 

f^W = (8 aS 8 rP , 8 ap 8 r5 ) , f uba ? = ((r>T b ) ap , (T b T') ap , tx[T a T b ] 8„ p ) . (3.19) 

Using Eq. ( ]3.17| ), we can rewrite Eq. (3A) in its final form 

P N 

-2tff = / W&Pi °+-.(.-) (P 1 ' • • • >Pn)C+~(—) (Pi,---,Pn), 
■' i=\ 

where the final operators are 



(3.20) 



0+ = ()■'■ "- fUi-aN . ( 3_2i) 
If we now use the same basis as in the usual color decomposition of the amplitudes 

MjCD (g a i (px)g a 2(p2) ■ -qTiPN)) = i^f k U ' [a2 - aN A k (l + ,2-,. ■ ■ ,N±) , (3.22) 



the MS Wilson coefficients are equal to the IR-finite parts of the color-ordered amplitudes to all 
orders 

C k + _. i .._ ) (p l ,p 2 ,. . . ,p N ) = A| n (l+,2-, . . . ,JV+) . (3.23) 



6 



Combining fixed-order helicity amplitudes with resummation using SCET 



Frank J. Tackmann 



To give a concrete and slightly nontrivial example, consider ggqqH. There are a total of six 
independent helicity operators, 

1 , 



a 



o 



o 



fibafi 
++(±) 
abap 
+-(±) 
abap 
-(±) 



'2+ • / 34± 



H 5 , 



jap 



/ 34± // 5 
r «j3 



(3.24) 



where the factors of 1/2 are symmetry factors. Using the above color basis for ggqq, the color 
decomposition of the QCD helicity amplitudes into partial amplitudes reads 



^(gl±g2±q3 + <j4-H 5 ) = i £ [T a °WT a °(2)] A((T(1 ± ),(T(2 ± );3+ 4J;5 H 



- i tr[T fll T" 2 ] 5 a3 « 4 fl( 1 ± , 2 ± ; 3+ , 4^ ; 5 ff ) . 



(3.25) 



The B amplitudes vanish at tree-level. From Eq. ( 3.25[ ) we can read off the Wilson coefficients, 

^(^,2^,3+ 4J;5») 



C±±f + )(pi,P2\P3,P4lP5) 



(3.26) 



A fin (2 ± ,l ± ;3+,4,r;5 // ) | . 
^(1+2^3+ 4r;5 H ) 

The Wilson coefficients for the negative helicity quark current can be obtained using charge conju- 
gation invariance, 

/ -1 

CxM-){puP2;P3,P4;Ps) =yCxM + ){puP2\PA,py,P5) where V = 



-10 | . (3.27) 
0-1 



4. Conclusions 



The tools to combine generic NLO calculations with a NNLL resummation are available. The 
required matching coefficients from matching QCD onto SCET contain the process-dependent hard 
virtual corrections, but do not depend on the choice of jet-resolution variable. We have shown how 
to explicitly construct an operator basis in SCET such that the IR-finite parts of the color-ordered 
partial amplitudes directly determine the hard machting coefficients. This provides a seamless 
interface to combine existing NLO calculations with higher-order logarithmic resummation. 

It is important to note that in general the different color structures in O f mix under renormal- 
ization. Except in the simplest cases, this already happens at NLL. Including the RGE running of 
the hard coefficients, the resummed cross section has the color structure 



a N 



H 



B a B b Hjj)®S 



U H -C, 



(4.1) 



where the hard evolution factors Uh as well as the soft function S are matrices in color space, which 
depend on the kinematics. In practice, this means that in order to perform the resummation one 
requires the NLO calculation to provide access to the individual color-ordered amplitudes. For this 
purpose it would be important to establish a common set of conventions for the helicity and color 
bases. 
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